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Abstract 

We provide the exact generating function for semi-flexible and super-flexible interacting 
partially directed walks and also analyse the solution in detail. We demonstrate that while 
fully flexible walks have a collapse transition that is second order and obeys tricritical scaling, 
once positive stiffness is introduced the collapse transition becomes first order. This confirms 
a recent conjecture based on numerical results. We note that the addition of an horizontal 
force in either case does not affect the order of the transition. In the opposite case where 
stiffness is discouraged by the energy potential introduced, which we denote the super-flexible 
case, the transition also changes, though more subtly, with the crossover exponent remaining 
unmoved from the neutral case but the entropic exponents changing. 

1 Introduction 

The collapse transition of an isolated polymer has continued to attract both theoretical and 
experimental attention. The canonical lattice model of single polymer collapse has been the self- 
avoiding walk with the addition of attractive potentials between non-bonded nearest-neighbour 
sites of the walk. This is known as the Interacting Self-avoiding Walk (IS AW). This model has 
yielded many important theoretical aspects of the physical problem though it is not exactly 
solved in the sense that the generating function of partition functions has not been explicitly 
calculated, in two or three dimensions. An exactly solved version of the model does exist however 
when the restriction of partial directness is imposed on the configurations of the self-avoiding 
walk in two dimensions. The model has been shown to display a tricritical-like collapse transition 
[I] as is predicted for the unrestricted model, though with different exponents. 

The Interacting Partially Directed Self-avoiding Walk (IPDSAW) model, and a closely related 
semi-continuous variant, on the square lattice was studied extensively in the early 1990's [21 El 
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HI El dj E] . It was noticed that this problem is in a family of related problems including lattice 
models of vesicles [TJ |8] whose solution can be written in terms of g-Bessel functions: moreover, 
direct correspondences occur between various models. Importantly, key work associated with 
the asymptotic analysis of the functions that arise in this class of problems was also completed 
[9]. Taken together these works completely solve and analyse the generating function, and free 
energy, of the IPDSAW model. In particular, the location of the collapse transition was found 
by Binder et al. [2] while the exact generating function was found by Brak et al. [3j in terms of 
g-Bessel functions. The tricritical nature of the collapse transition was elucidated by Owczarek 
et al. [1] and the full asymptotics of the generating function can be deduced from the work of 
Prellberg [9]. 

The addition of a stiffness parameter to mimic the effects of persistence length [10] and a 
stretching parameter to model the pulling of a polymer by an external force has more recently 
been studied in the context of the ISAW model [TT]. While a parameter called a pulling force 
was not explicitly mentioned in the work on the IPDSAW it was implicitly part of the set up of 
the model, as we shall see below, since the horizontal and vertical steps of the walk were given 
separate fugacities. It was shown [HE] that differentiating the horizontal and vertical fugacities 
does not affect the nature of the collapse transition. Separate analysis of the IPDSAW model 
with the force interpretation being explicit confirms this [12] . On the other hand, the addition 
of a stiffness parameter so that the polymer is now semi-flexible was not included in the original 
definition of the model. The IPDSAW has recently been reconsidered by Zhou et al. [13| . 
Interestingly, they conjectured, on the basis of Monte Carlo simulation, and an approximation 
scheme allowing precise numerical estimates of thermodynamic quantities, of semi-stiff IPDSAW, 
that positive stiffness changes the order of the collapse transition to first-order. The three- 
dimensional semi-stiff ISAW model has been studied by Grassberger and Hegger |10] some time 
ago and they showed that the collapse transition does indeed become first order though only 
for a finite amount of applied stiffness. That is, small stiffness parameter values do not change 
the nature of the collapse transition. A related model shows similar behaviour [T3]. In this 
paper we solve exactly the IPDSAW with stiffness parameter, which we shall now refer to as the 
Variably-Flexible Interacting Partially Directed Walk (VFIPDSAW) and analyse the model in 
the full parameter space. We show that not only does the collapse transition become first order 
when the stiffness parameter is positive (semi-flexible case) but it is also modified, though still 
tricritical, when the stiffness parameter is negative (super-flexible case). 

1.1 The model 

Consider the square lattice and a self-avoiding walk that has one end fixed at the origin on that 
lattice. Now restrict the configurations considered to self-avoiding walks such that starting at 
the origin only steps in the (1,0), (0, 1) and (0, —1) are permitted: such a walk is known as a 
Partially Directed Self-avoiding Walk (PDSAW). For convenience, we consider walks that have 
at least one horizontal step. Let the total number of steps in the walk be L and the number of 
horizontal steps be N. Hence, we have L > N > 1. An example configuration along with the 
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associated variables of our model are illustrated in Figured! We begin by recalling the definition 
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Figure 1: An example of a Partially Directed Walk of length L = 21 and horizontal length 
N = 8 (the bold black path) with the four parameters x associated with horizontal steps, y 
associated with vertical steps, lo associated with nearest-neighbour 'contacts' (shown as (red) 
intertwined curves) and the stiffness parameter a associated with the sites (highlighted in grey 
(green)) between two consecutive horizontal steps . The weight of the configuration shown is 
x 8 y 13 oj 6 a 2 . 

of the IPDSAW model and then add the "stiffness parameter" , the addition of which defines 
our model. 

IPDSAW To define the IPDSAW model we add various energies to properties of this walk 
and hence Boltzmann weights to the walk. Firstly, any two occupied sites of the walk not 
adjacent in the walk though adjacent on the lattice are denoted "nearest-neighbour contacts": 
see Figure [U An energy — J is added for each such contact. We define a Boltzmann weight 
u) = e l3J associated with these contacts, where j3 = l/fegT and T is the absolute temperature. 
An external horizontal force / pulling at the other end of the walk adds a Boltzmann weight p N 
and p = e@f a ° , with ao being the length of a lattice bond. The partition function Z| PDSAW (o;,p) 
of the IPDSAW model is 

zr sAW (^ P ) = £ u/yv, (i.i) 

PDSAW of length L 

where m is the number of nearest-neighbour contacts in the PDSAW. The generating function 
G ipdsaw (z,lo, P ) 

oo 

G IPDSAW (z, W>P ) = £ Zl PDSAW (u,p)z L , (1.2) 

L=l 

so z can be considered as fugacity for the steps of the walk and the generating function a 
"generalised partition function" [1J. In previous work [I] an alternate generating function 
G IPDSAW {x,y,u) was considered, where instead of a force parameter p horizontal steps were 
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weighted with a fugacity x while vertical steps were given a fugacity y: see Figured) Hence we 
have 

G IPDSAW (z, W ,p) = G IPDSAW (pz,z,u,) . (1.3) 

Clearly considering a separate horizontal fugacity is equivalent to considering a horizontal pulling 
force at the level of generating functions. 

VFIPDSAW To define the VFIPDSAW we now add an energy —A to each site between 
consecutive horizontal steps of the walk: see figure [TJ Note that for A > consecutive hori- 
zontal steps are favoured and so this is the positive stiffness, or semi-flexible, regime while for 
A < consecutive horizontal steps are discouraged so this is the negative stiffness, or super- 
flexible regime. If £ is the number of such "stiffness" sites in a particular IPDSAW then such 
a configuration is associated with an additional Boltzmann factor a e where a = e^. That is, 
each configuration has weight u m p cr. Note that one could have equivalently chosen to weight 
every bend, or change of direction of the walk, with a weight b say. That is, each configuration 
has weight u> m p N b k . However, since the number of such bends k is related to the number of 
horizontal straight segments, £, assuming for convenience at least one vertical step in the walk, 

, = <?£z*z!), (1 . 4) 

then substituting a = 1/b 2 and setting p to pb 2 gives the same weight for each configuration 
(barring an overall factor of b). 

The partition function for VFIPDSAW is defined as 

z VFIPDSAW (u;iP;(T)= £ wVVi (L5) 

PDSAW of length L 

while the generating function, analogously to the fully-flexible case above, C VFIPDSAW (x, y,u,a), 
is given as 

oo 

G VFlPDSAW (pz,z,u;,a) = £ Z™ SAW \u;,p,a)z L . (1.6) 

L=l 

Clearly one can recover the fully flexible case by setting a = 1: 

G 1PDSAW (pz,z,u J ) = G VFlPDSAW (pz,z,u J ,l). (1.7) 

Setup As we consider the VFIPDSAW model from now on we shall drop the superscript 
VFIPDSAW. The singularity structure of the generating function as function of z determines 
the free energy. The reduced free energy is defined as 

n(u),p,o-) = - lim -j-log [Z L (uj,p,a)] (1.8) 

L— >oo Li 

and is given by 

k(uj, p, a) = log z s (u,p, a) , (1.9) 
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where z s (uj,p,a) is the closest singularity (on the positive real axis) of the generating function 
G(pz, z, uj, a) in the variable z to the origin. Note also that 



Z L (u,p,a) = [z L ]G(pz,z,uj,a) = j> G(pz,z,u),o) 



dz 

yL+l ■ 



(1.10) 



In order to find the generating function it is advantageous to rewrite it in the following 
way. We can describe the PDSAW configurations in a natural way through the length rj of 
vertical segments between two horizontal steps, measured in the positive y-direction. Each 
PDSAW begins with a vertical segment of height r\ followed by an horizontal step. Thus, we 
associate to each configuration an ./V-tuple (ri,r2, . . . , rjv) corresponding to a configuration of 
total length L = J2iLi \ r i\ + N. The energy due to the nearest-neighbour contacts for each of 
these configurations is then 

- J u(n,r 2 , ...,r N ), 



where 



N-l 



u(ri,r 2 , ...,r N ) = min (N> |n+i|)©(-^^+i) , 



(1.11) 
(1.12) 



i=i 



where 0(r) is the Heaviside step function: 

G(r) = { 



r < 0, 
1/2 r = 0, 

1 r > 0. 



(1.13) 



The number of "stiffness sites" £ is then given by the number of times r i = for any 1 < i < N. 
We get the generating function by summing over all possible lengths as 



oo L 



G(x,y,u,a) = J^JZ 



N 



E 



y L-N UJ u(r 1 ,r2,...,r N ) .e ^ 



that is, 



L=1N=1 \ ri \ + \ r2 \ + ... + \ rN \ =L - N 



G(x,y,u;,a) = Y,x N Y.y M E 

JV=1 M=0 |ri|+|r 2 |+...+|rjv|=M 



(1.14) 



(1.15) 



2 Exact solution of the generating function 

In order to derive an expression for G(x,y,uj,a), consider the generalised partition functions 
G r = G r (x, y, uj, a) for walks that start with a vertical segment of height r, so that 



G(x,y,u,a) = ^G r . 
Then we can concatenate these walks to get a recursion relation for G r as follows: 

G r - - xy |r| { 1 + a 5r >°G + ^ uj u ^ s) G s 

s&I\{0} 



(2.1) 



(2.2) 
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It follows that 

G = x{l + (a-l)G + G} , (2.3) 

so that 

G = xu{l + G}, (2.4) 

where 

u=- } -. (2.5) 

1 — x(a — 1) 

Using the symmetry G r = G- r and then restricting to r > 0, we can further simplify to 

G r = xy r jl + a 5 ^G + f^l + oo min ^)G s ^ , (2.6) 

which will be the starting point of our investigation. Since r = is now special we will need to 
consider r = 1 separately also: 



Gi =xy|l + G + (l+cu)^G s \ . (2.7) 

Now using 



s=l 



^2G S = -(G-G ) (2i 



s=l 

gives 



G^xyU + ^G^^G}. (2.9) 



Now using equation (|2.4|) we obtain 



Hence the ratio G\/Gq can be written in terms of 1 + G. By solving for (1 + G) one finds 

(2.11) 



_yG 

We will now derive a homogeneous second order difference equation which we can solve 
using the same ansatz used previously |l!J. Using the scaling behaviour of the solutions, we can 
eliminate one of the two linearly independent solutions. We then write the general solution of 
(12.6H as an expression involving the quotient of two g-hypergeometric functions. 

Taking differences in (|2.6[) . we first eliminate the inhomogeneous term, 

_. oo 

G r+1 - yG r = 5 rfi xy(l - a)G + xq r+1 (l - -) V G s . (2.12) 

s=r+l 

Here, we introduced for convenience the new variable q = yuj. Upon taking differences a second 
time, we are left with 

(G r+2 - yG r+l ) - q(G r+1 - yG r ) = -5 rm xy{\ - a)G - xq r+2 (l - -)G r+l . (2.13) 
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We now solve this equation for r > 1 and subsequently solve for r = 0. In the case of no 
interaction (u = 1), the right hand side of this equation is zero (for r > 1) and we have a simple 
homogeneous difference equation with constant coefficients. Its characteristic polynomial P(l) 
is 

P(l) = (l-y)(l- q ) (2.14) 

and the solution is given by G r = A\y T + A 2 q r ■ 
This motivates the ansatz |15j 

oo 

G r = r^2q nr c n , (2.15) 

n=0 

with c n = c n (x,g,ti;) independent of r > 1, which inserted into (|2.13|) gives 

P(/)c + f; g nr C P(/cf )c„ + xq{\ - ^)k n c n -i) = 0. (2.16) 

n=l ^ ' 

This equation is solved by 

P(l) = 0, i.e. h = y and l 2 = q, (2.17) 

and, choosing cq = 1, 

" -*<?(! - _ (-xw(l - £)Z)VS) 



Here we have used the standard notation 

n 

(x;q)n= Hil-xq™- 1 ). (2.19) 



m=l 

Defining 



fi)(-ty 



H{y,q,t) = Y j ( q A \ , (2-20) 
we now can write the general solution of (I2.13L for r > 1, as 

G r = A lV r H(y, q, x(l - -)q 1+r ) + A 2 q r H(qu, q, xu(l - -)g 1+r ). (2.21) 
We remark that the function H is directly related to a basic hypergeometric function [16] 

H(y,q,t) = 1 M0,y,q,t), (2-22) 

which can be seen to be a limiting function of 2^1 and that is the ^-deformation of the more 
familiar hypergeometric function 2 Fi. Analogously, the function H can be understood (apart 
from some normalising factors and seen by taking the limit q —* 1) as a g-generalisation of Bessel 
functions. One can easily verify that H(y, q, t) satisfies the following recurrence 

qH(y, q, t) + (tq - (y + q))H(y, q, qt) + yH(y, q, q 2 t) = . (2.23) 
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Returning to the analysis we see that, for \q\ < 1, H(y,q,tq r ) is uniformly bounded in r, so 
that we can write 



\G r \ < const (q r + y r ). 
This we insert into (|2.6p and, assuming < uj 2 y < 1 < u we get 

/ i — 1 <y 



(2.24) 



\G r \ < const y r 1 + ^(ujq) 3 + uf ^ <f 



< const y r (l + (wg) r ) < const y r . 



(2.25) 
(2.26) 



As -ff (y, g, tg ,r ) — > 1 for r — » oo, we see that in fact Ai = 0. The reason for this is that we obtained 
the homogeneous difference equation (|2,13p by taking differences from (j2.6j) . thus introducing 
additional solutions. 

We now note that the ratios GV+i/GV contain no unknown constants. In fact, defining 

H(y,q,qt) 



H(y,q,t) ' 



we find 



^=H(y,q,x(l--) q 1+r ) forr>l. 
yG r <jj 



(2.27) 



(2.28) 



Note in passing that successive ratios are related via the following recurrence for 7i{y, q, t) derived 
from equation (I2.23|) . 

H(y, q,t)=q[y + q-tq- yH(y, q, qt)}' 1 . (2.29) 

In fact, the ratio G\/Gq is given by a very similar expression. For r = 0, the recursion ()2. 131) 
can be rewritten as 

(G 2 - yGr) - q{G 1 - ^G ) = -xq 2 (l - -)G l , (2.30) 

U UJ 

from whence one can conclude that 



= W(y,g,x(l - -)q) . 
yG uj 

Inserting this ratio into equation (|2.1ip and substituting q = yu>, we have 



(2.31) 



l + G 



(1-w) 



H(y,yuj,x(uj - l)y) 



1 xu 



-1 



(2.32) 



We note immediately that the stiffness parameter a only enters (via u) in one term of this 
expression. For a = 1 (i.e. u = 1), this is exactly equation (4.27) in pQ. 

Our final expression for the solution of the generating function for the Variably Flexible 
Interacting Partially Directed Walk in the variables z,uj,p,a is therefore 

-i 



l + G 



(1-u,) 



Ti(z, zu:,pz (lu — 1)) 



1 + | (l-o>) 



pz 



1 — pz(a — 1) 



(2.33) 
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As we will see below the case q = 1 is important. From (|2.29p it follows that H(y, 1, t) is the 
root of a quadratic equation, and so 

1 



H(y,l,t) 



2y L 



l + y-t-y/(l + y-t) 



4y 



(2.34) 



where the branch has been chosen such that Tl(y,l,t) = 1 + 0(t). There is an algebraic 
singularity at (1 + t — y) 2 = 4y and the solution is real-valued as long as (1 + 1 — y) 2 > 4y. 
One can therefore solve for the generating function along the curve yuo = 1 as 



1 + G 



(1-u) 



n(-,i,x(i--)) 

00 U) 



(i + w) + 0_ 



-xu 



(2.35) 



which is now algebraic. 



3 Analysis of the Phase diagram 
3.1 General considerations 

One can immediately observe that the generating function G has singularities at the singularities 
of Tt(y,yu!,x(oj — l)y) and when the denominator is zero, that is, at solutions u> c (x,y,u) of 

H(y, yu, x(uj - l)y) = + . (3.1) 

There is an essential singularity of TC(y,yuj,x(l — oS)y) when yco = 1. On the other hand, when 
the denominator is zero, G has a pole, and the locus u c (x, y, u) of this pole depends analytically 
on the parameters as long as yu < 1. If there is no zero of the denominator for yuj < 1, then the 
closest singularity is given by the essential singularity of H(y, yu>, x(l — oj)y) at yu = 1 where the 
generating function converges. On this curve, we obtain from (12.350 that there is an algebraic 
singularity at 



W p (x,MJ = -^j . (-J--JJ 



and for u > 1 a simple pole at 

(1 + -ux)(l + 2x — ux) 
(1 — ux)(l — 2x + ux) 

These singularities coincide when u = 1, at which value the nature of the algebraic singularity 
changes. Note that for u < 1 the pole disappears. 

As stated above for any fixed to the generating function as a function of y either has a pole 
given by the solution of u = uj c (x,y,u) or has a singularity on the curve y = 1/lo. Therefore, 
for u < 1, uj c (x,y,u) meets the algebraic singularity u a (x) as the curve yu = 1 is approached, 
whereas for u > 1, oj c (x,y,u) meets the pole uj p (x,u) as the curve yu = 1 is approached. To 
investigate the behaviour of G near the curve yu) = 1 more closely, we need the asymptotic 
behaviour of Tt(y, q, t) as q — ► 1. Away from the algebraic singularity, i.e. for (1 + i — y) 2 > 4y, 
we can use (|2.29p to derive an asymptotic expansion in e = 1 — q, 

oo 

H(y,l-e,t)^J2n (n) (y,t)e n (3.4) 

?1=0 
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with the first terms given by TC^(y, t) = 7i(y, 1, t) and 



yn(°)(y,tf 



1 



1 + ytH^ ( y , t) (y,t) + (t- (y , i) 



1 I | - y+t ~ l 
2y I ^(l+y-t)2_4y 



+ 



2yt 



l+y-t) 2 -4y 



(3.5) 
(3.6) 



Close to the algebraic singularity at q = 1, the singularity structure is significantly more compli- 
cated, but has been thoroughly elucidated in [9j. Using Lemma 4.3 from [S], a result completely 
analogous to Theorem 5.3 in [9] can be obtained for H(x,q,t), i.e. an asymptotic expansion in 
q = 1 — e uniformly valid for all values of t and x, which reads 



H{y,l-e,t) 



1 

2^ 



1 + y-t 



Ai'(ae- 2 / 3 ) 



a l/2 e -l/3 Ai(ae" 2 / 3 ) 



^(l + y-ty-Ay 



:i + 0(e)). 



(3.7) 

Note that for e — > the expression multiplying the square root in (13. 7p tends to 1 as is necessary. 
Here, a = a(y,t) is a function of y and i which is known exactly [S]. While the precise form of 
a is rather cumbersome, it simplifies considerably near the critical point, and we find 

4/3 {1 + y _ t) 2_ Ay 



a(y,t) 



i-{t-yf 



(3.8) 



for small (1 + y — t) 2 — 4y. This implies that here 



W(y,l-e,t) ~ — 
2y 



1/3 Ai'(q e - 2 / 3 )(l-(t-y) 2 ) 2 / 3 
+ 2/ + Ai(ae" 2 / 3 ) 2V3 



(3.9) 



The behaviour of this expression is determined by the function /(#) = — Ai'(z) / Ai(z), the graph 
of which is shown in figure [21 The large- z asymptotics allows for matching for e — * and positive 
a. For negative a, the argument of / is negative. As f(z) has a simple pole at z = —2.3381 . . ., 
for any fixed a < we have a pole at a finite value of e. As a tends to zero, the locus of this 
pole scales as e 2 / 3 . 



3.2 Fully flexible case (A = 0) 



This case is the one considered in our earlier work [lj. Here u = 1, and the simple pole at 
u c (x, y, 1) approaches an algebraic singularity at uo t = ( x=f ) > 1, at which the generating 



function diverges. In particular, we find 

Ai'(ae- 2 / 3 ) 



-i 



1 + G 



LO - 1 



n 



1 /2 e -i/3 Ai(Qe -2/3) J ^(l +w + x(l-o;)) 2 -4 W 



(3.10) 



as e = 1 — yw — > 0, with a = a(— , x(l — — )) given by f|3.8|) . Near the transition, a is small and 
we can write 



+ Ai'(ae- 2 / 3 ) 



(3.11) 
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Figure 2: The graph of f{z) = — Ai'(z)/ Ai(z). The function f has a simple pole at z = 
—2.3381 . . ., a zero at z = —1.0187. . ., and is asymptotic to ^fz for large z, which is plotted for 
comparison. 



with A — {u — 1)/V 2 - (x(uj - 1) - l) 2 ] 2 / 3 (2o;) 1 / 3 . For u = u t , a = and G diverges as e -1 / 3 . 
For uj > ut, a > and G tends to a finite value given by (|2.35p . For to < u>t, a < and G 
has a simple pole at ae 2 / 3 ~ —1.0187. . .. We accordingly have a second-order phase transition 
characterised by 

11,2 . 
7«=2 7t = 3 0=3, (3-12) 

with the exponents as defined by Owczarek et al. pp. 

Changing to the variables z, tv,p, a this can be formulated as follows. This is the case a = 1, 
and there is a curve of simple poles given by uj c (pz,z, 1) approaching the curve loz = 1 at Ut 
given by u>t = ( jjjz^ J ■ For p = 1, the solution is u>t = 3.3829. . .. Near the transition, (|3.10p 
holds with the appropriate substitution x = pz and y = z. While the location uj c (pz, z,l) of the 
poles, as well as ojt change as a function of p, the character of the phase transition does not. 

3.3 Super-flexible case (A < 0) 

Now u < 1, and the simple pole at u> c (x,y,u) approaches an algebraic singularity at oJt = 
(l^x) > 1) a t which the generating function converges. In particular, we find 

1 + G 7 a-// \ ^^^^=^^^^= ( 3 -l 3 ) 

x(u - 1)(1 - u) + (- n ^i7Jj-yJ Va+^ + ^l-^)) 2 -^ 



as e = 1 — yuj — > 0, with a = a(— , x(l — ^)) given by (|3.8|) . Near the transition, a is small and 
we can write 

1 + G _ 2/3 , (3.14) 
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with A as above. For u = Ut, a = and G converges with the singular part scaling as e 1//3 . For 
u > cot, a > and G tends to a finite value given by (|2.35p . For co < u>t, a < and G has a 
simple pole at some value of e where - Ai'(z)/ Ai(z) < 0, i.e. -2.3381 .. . < ae 2 / 3 < -1.0187.... 
We now find a second-order phase transition with 

1 1,2 . 

7u = - g 7t = - g ^ = 3 • ( - 3 - 15 ' ) 

As u — > 1, we recover the fully flexible case discussed above. 

Changing to the variables z, uj,p, a this can be formulated as follows. This is the case a < 1, 
and there is a curve of simple poles given by uJ c (pz, z, (1 — pz(a — 1)) _1 ) approaching the curve 

2 



ujz = 1 at cot given by uJt — ( ) • This transition point therefore is independent of the value 



K p-u> tj 

of a, and is identical to the one obtained in the fully flexible case. Near the transition, (|3.13p 
holds with the appropriate substitution x = pz, y = z, and u = 1/(1 — pz(a — 1)). While the 
location u> c (pz,z, (1 — pz(a — l)) -1 ) of the poles, as well as u>t change as a function of p, the 
character of the phase transition does not. 



3.4 Semi-flexible case (A > 0) 

(1— ux)(l— 2x+ux) 



Now u > 1, and the simple pole at lo c (x, y, u) approaches a simple pole at uj t = i!^"!!!^ 1 , ux \ > 



1. In particular, we find near the transition that 

l + G = = 1 — : — (3.16) 

x(lo - 1)(1 - u) + +co + x(l -co)) 2 -4u - 2eW( 1 )(i,x(w - 1)1) 



ase = l-yu -> 0, with^^ 1 ) given by (|33]) . Note that x{uj-l){l-u)+y/{l + u + x{l - uo)) 2 - 4to 
is asymptotically linear in u — Ut and is negative for u < uot- Note further that < 0. For 
lu > u>t, G tends to a finite value as e — > 0. This expression diverges with a simple pole as u 
approaches uJt- Similarly, for oj = w t , G diverges as e -1 . For uj < u> t , G has a simple pole at 



_ s(u;-1)(1-m) + ^(l + cv + xjl-u)) 2 -^ 

2H« ( I^-1)1) • (3 - 17) 

We accordingly find a first-order phase transition with 

7u = l lt = l = 1. (3.18) 

As it — ► 1, T^ 1 ' diverges, and the character of the phase transition changes as characterised 
above. 

Changing to the variables z, u,p, a this can be formulated as follows. This is the case a > 1, 
and there is a curve of simple poles given by u) c (pz, z, (1 — pz(a — 1)) _1 ) approaching the curve 
ujz = 1 at ujt given by 

= (2(1 - a)p 2 + (2 - a)uj t p + uj 2 ){(2 - a)p + uj t ) ^ 
(2(cr - l)p 2 - w t ptr + ujf)(uj t - per) 
While this is quite cumbersome in general, some special values have simple solutions. For 
example, at p = 1 and a = 2 we find uj t = 2 + V2 = 3.4142 .... Near the transition, (^T6|) holds 
with the appropriate substitution x = pz, y = z, and u = 1/(1 — pz[a — 1)). While the location 
uj c (pz, z, (1 — pzio — 1)) _1 ) of the poles, as well as uj t change as a function of p, the character of 
the phase transition does not. 



12 



4 Conclusion 

We have analysed the exact solution of a two-dimensional lattice model of a single polymer in 
solution containing parameters that vary the intra-polymer attraction, the amount of horizontal 
stretching force applied and the amount of stiffness. The restriction of partial directness is 
required to ensure solvability. We find that a tricritical collapse transition takes place for no 
stiffness or negative stiffness, and that this is unaffected by an horizontal force. The entropic 
exponents are different in the negative stiffness regime to those in the zero stiffness regime. 
On the other hand, when the polymer becomes semi-stiff the collapse transition immediately 
becomes first-order. 
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